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We calculate entropy-temperature curves for interacting bosons in unit filled optical lattices for both homo- 
geneous and harmonically trapped situations, and use them to understand how adiabatic changes in the lattice 
depth affect the temperature of the system. In a translationally invariant lattice, the zero tunneling limit facili- 
tates a rather detailed analytic description. Unlike the non-interacting bosonic system which is always cooled 
upon adiabatic loading for low enough initial temperature, the change in the excitation spectrum induced by 
interactions can lead to heating. Finite tunneling helps to reduce this heating. Finally, we study the spatially 
inhomogeneous system confined in a parabolic potential and show that the presence of the trap can significantly 
reduce the final available temperature, due to the non- vanishing superfluid component at the edge of the cloud 
which is present in trapped systems. 



I. INTRODUCTION 

Cold atoms in optical lattices provide a system for realiz- 
ing interacting many-body systems in essentially defect free 
lattices jjl, and have been an active area of research in re- 
cent years. The strong interest in this system is due in part 
to the ability to dynamically control lattice parameters at a 
level unavailable in more traditional condensed matter sys- 
tems. Lattice-based systems are typically governed by three 
sets of energy scales: interaction energies U, tunneling rates ,/ 
and the temperature T. In atomic systems, the energies U and 
J can be controlled by adjusting the lattice, and their values 
can be measured and/or calculated easily. Unlike condensed 
matter systems, however, it is experimentally difficult to mea- 
sure very low temperatures, (kT < J, kT < U, here k is the 
Boltzmann constant), and the temperature has so far only been 
inferred in a few casesj^-IHBSSl- Absent good thermome- 
ters, and given the ability to dynamically change the density 
of states, it is important to understand the thermodynamics of 
experimentally realistic systems in order to estimate the tem- 
perature. 

It has been pointed out that loading sufficiently cold, non- 
interacting atoms into an optical lattice can lead to adiabatic 
cooling, but the cooling available in a real system will 

clearly depend on and be limited by interactions. It can also 
depend on the (typically harmonic) trapping potential, which 
provides an additional energy in the problem, as well as on 
the finite size of the sample. Here, we calculate the entropy 
of bosons in unit filled optical lattices for homogeneous and 
trapped cases. We provide good approximate, analytical ex- 
pressions for the entropy for various cases, including finite 
number effects which allow for comparison of temperatures 
for adiabatic changes in the lattice. 

For translationally invariant lattices at commensurate fill- 
ing, the reduced density of states associated with the gap that 



appears in the insulating state presents a significant limitation 
to the final temperature when raising the lattice f3\. The pres- 
ence of the trap, and the associated superfluid-like component 
at the edges can significantly increase the density of states, 
however, allowing for lower final temperatures. 

In this paper we make the assumption of adiabatic load- 
ing and thus calculate the lowest possible final temperature 
achievable from a given initial temperature during the loading 
process. We realize that to be fully adiabatic might be experi- 
mentally challenging, however our calculations could be used 
to benchmark the effect of the loading on the temperature of 
the atomic sample. 

The paper is organized as follows: We start by introduc- 
ing the model Hamiltonian and our notation. In Sec. Illllwe 
focus on the translationally invariant case. We first develop 
analytic expression for the thermodynamic quantities in the 
J = limit and then we use them to calculate the final tem- 
perature of the atomic sample assuming we start with a dilute 
weakly interacting BEC, described using the Bogoliubov ap- 
proximation. Next we study how finite size effects and finite 
J corrections modify the final temperature of the sample. In 
Sec. llVl we discuss the effects of a spatial inhomogeneity in- 
duced by an additional parabolic potential and finally in Sec. 
|V]we conclude. 

II. BOSE-HUBBARD HAMILTONIAN 

The Bose-Hubbard {BH) Hamiltonian describes interacting 
bosons in a periodic lattice potential when the lattice is loaded 
such that only the lowest vibrational level of each lattice 
site is occupied and tunneling occurs only between nearest- 
neighbors Ql] 

(iJ) j 
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Here aj is the bosonic annihilation operator of a particle at site 
j = {jx, jy, jz}, fij = ajoj, and the sum (i, j) is over nearest 
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neighbor sites. U is the interaction energy cost for having two 
atoms at the same lattice site which is proportional to the scat- 
tering length as, Vj accounts for any other external potential 
such as the parabolic magnetic confinement present in most 
of the experiments and Jj j is the hopping matrix element be- 
tween nearest neighboring lattice sites. 

For sinusoidal separable lattice potentials with depths 
{Vx,Vy,Vz} in the different directions, the nearest neigh- 
bors hopping matrix elements, {Jx,Jy, Jz}, decrease expo- 
nentially with the lattice depth in the respective direction and 
U increases as a power law; U oc asiV^VyVzY^'^^. 

III. HOMOGENEOUS LATTICE 
A. Thermodynamic properties in the J = limit 

In this section we calculate expressions for the thermody- 
namic properties of strongly correlated bosons in a spa- 
tially homogeneous lattice (Vi = 0), with M sites. For the 
case where Jx.y.z = 0, (relevant for very deep lattices) the en- 
tropy can be calculated from a straightforward accounting of 
occupation of Fock states, and is independent of the number 
of spatial dimensions. We derive expressions for the entropy 
per particle as a function of AI, N, U and the temperature T, 
in the thermodynamic limit where N ^ oo and M ^ oo, 
while the filling factor N/M remains constant. 

In the Jx,y,z = limit, Fock number states are eigenstates 
of the Hamiltonian and the partition function Z can be written 



Z{N,M) = n{nr)e- 



(2) 



where (3 = (kT) ^, k is the Boltzmann constant. Here we 
use the following notation: 

• The quantum numbers give the number of wells with 
r atoms, r = 0, 1, . . . , iV, in a particular Fock state 
of the system. For example for a unit filled lattice the 
state |1, 1, . . . , 1, 1) has quantum numbers rii = N and 

rir^ii = 0. 

. Er = ^r{r - 1) 

• The sum is over all different configurations {rir} which 
satisfy the constrains of having N atoms in M wells; 



A/, 



rrir = N, 



r=0 

E 

r=0 



(3) 
(4) 



• ri(nr) accounts for the number of Fock states which 
have the same quantum numbers and thus are degen- 
erate due to the translational invariance of the system: 



Notice that without the particle number constraint, Eq.0 and 
Eq. ^ could be easily evaluated. It would just be given by 



J2ziM,N) 



N 
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no,"!,--- 



Ml 



/3-Bl^nl 



M 



-I3E,. 



(6) 



However, the constraint of having exactly N atoms, Eq.Q, in- 
troduces some complication. To evaluate the constrained sum 
we follow the standard procedure and go from a Canonical to 
a grandcanonical formulation of the problem. 
Defining the grandcanonical partition function: 



M 



S{M) = Y,Z{N',M)e 



N' 



(7) 

and using the fact that S{M) is a very sharply peaked func- 
tion, the sum in Eq.Qi can be evaluated as the maximum value 
of the summand multiplied by a width AA^*: 

E:{M) ^ Z(N,M)e'^^'^AN*, (8) 

Taking the logarithm of the above equation and neglecting 
the term ln(AAf*), which in the thermodynamic limit is very 
small compared to the others (AA^* <C N), one gets an 
excellent approximation for the desired partition function, 

Z{N' = TV, A/): 



\n[Z{N,M)] = -iSi^iN + ln[S{M)]. 



(9) 



The parameter /i has to be chosen to maximize 
Z{N', Af )e'^''^' at N. This leads to the constraint: 



9 = 



(10) 



(11) 



where g — N/M is the filling factor of the lattice, Ur is the 
mean density of lattice sites with r atoms, and ^ is the chemi- 
cal potential of the gas. 

From Eg. ( II 1> and Eq.(0 one can calculate all the thermo- 
dynamic properties of the system. In particular, the entropy 
per particle of the system can be expressed as: 



S{M, N) = fc(-/3/i + 1 ln[^(Af)] + l3E), (12) 
where E ^ 1 /N Ern,. is the mean energy per particle. 



1. Unit filled lattice M ^ N 

For the case M = N li is possible to show that, to an ex- 
cellent approximation, the solution of Eg. (II H is given by: 



M 

no, ni, . . . ,n7v 



no!ni! . 



(5) 



(13) 
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with C = 1.432. Using this value of /i in the grandcanoni- 
cal partition function one can evaluate all the thermodynamic 
quantities. 

• Low temperature limit {kT < U) 

In the low temperature regime ^ ~ C//2. By replacing fjL = 
[7/2 in Eq.(0 one can write an analytic expression for S and 
E (and thus for 5 ) in terms of Elliptic Theta functions 

^?3(^,g) =l + 2Er=i'7"'cos[27iz]: 



E{N) 



pPU/2 



1 + 



E = 



U 



2 + (0,6-^^/2) 
2 + eW2[l + i?3(0,e-/3^/2)] 



(14) 



(15) 



with {z,q) = ddj, {z,q)/ dq. In this low temperature regime 
one can also write an analytic expression for 



2p-PU/2(r-lf 



2 + e'3^/2[l + ,93(0, 6-^*^/2)] 



(16) 



High temperature limit (kT > U) 



In the high temperature regime (3p ~ — ln[(l + g) / g] which 
is just /3yU = — In 2 for the unit filled case. This can be easily 
checked by setting /3 = in Eg. (II 1> and solving for /i. 

For large temperatures, /3 ^ , the grandcanonical parti- 
tion function and the energy approach an asymptotic value; 
\n[E{M)] M[\n{l + g)], E Ug. Therefore the 
entropy per particle reaches an asymptotic plateau S/k 



ln[a„] 



. This plateau can be understood 
is the number of all the possible 



N 

because Qo = ^M^-iy.m 
accessible states to the system in the one-band approxima- 
tion (total number of distinct ways to place bosons in M 
wells). It is important to emphasize however, that the one- 
band approximation is only valid for kT ^ Egap, where 
Egap is the energy gap to the second band. For example, 
for the case of '^^Rb atoms trapped in a cubic lattice potential. 



14 = K 



Vz, Egap > IOC/ for lattice depths 14 > '^Er. 



Here, Ej^ is the recoil energy, and E^ = h? / {%m<P) where 
d is the lattice constant and m the atoms' mass. At higher 
temperatures the second band starts to become populated and 
thus the model breaks down. 

In Fig n we plot the entropy per particle as a function of 
temperature for a unit filled lattice. The (red) dash-dotted line 
corresponds to the numerical solution of Eq. and Eq. ( II U . 
The solid line (barely distinguishable from the numerical so- 
lution) corresponds to entropy calculated using the analytic 
expression of p given in Eq. ( I13> . The (blue) dashed line cor- 
responds to the analytic expression of the entropy derived for 
the low temperature regime in terms of Elliptic Theta func- 
tions: Eqs.( I14> and dlSl l. From the plots one can see that Eq. 
( I13> is a very good approximation for the chemical potential. 
Also the analytic expression derived for the low temperature 




FIG. 1 : (color online) Top: Entropy per particle as a function of the 
temperature T (in units of U) for a unit filled lattice in the J = 
limit. Dash-dotted (red) line: Eq. <12> calculated using the numerical 
solution of Eqs. and <1U : Solid (black) line: entropy calculated 
using Eq. <13> for the chemical potential; Dashed (blue) line: Ea. <12> 
calculated using the low-temperature analytic solutions: Eqs. <14> 
and J15> . Bottom: Average occupation number Hr as a function of T 
(in units of U) . The conventions used are: 7T1 (continuous line), no ( 
dashed), n2 (dotted-dashed) , na ( crosses) and Hi (dots). 



regime reproduces well the numerical solution for tempera- 
tures kT < U. 

It is also interesting to note the plateau in the entropy ob- 
served at extremely low temperatures, kT < 0.05C/. This 
plateau is induced by the gapped excitation spectrum char- 
acteristic of an insulator which exponentially suppresses the 
population of excited states at very low temperatures. As we 
will discuss below the range of temperature over which the 
plateau exists is reduced if J is taken into account. 

In the Fig.l we also show n^, the average densities of sites 
with r atoms vs temperature calculated using Eqs.( I13> and 
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(II 1> . In particular ni is important because lattice based quan- 
tum information proposals 11(111 llfT^l rely on having exactly 
one atom per site to iniziaUze the quantum register and popu- 
lation of states with r 7^ 1 degrades the fidelity. Specifically 
we plot ni (solid line), tTq (dashed line), n2 (dotted-dashed), 
ris (crosses) and 7X4 (points). 

In the entropy-plateau region of Fig.l, corresponding to 
kT < 0.05U, particle-hole excitations are exponentially in- 
hibited and thus ni is almost one. For temperatures kT < 
U/2, fiQ is almost equal to n2, meaning that only particle- 
hole excitations are important. As the temperature increases, 
kT > U/2, states with three atoms per well start to become 
populated and therefore rTg becomes greater than 712 ■ The 
population of states with r > 3 explains the break down of 
the analytic solution written in terms of elliptic functions for 
kT > U/2 as this solution assumes no —ri2. For kT > 2U, 
even states with 4 atoms per well become populated and the 
fidelity of having unit filled wells degrades to less than 60%. 

B. Adiabatic Loading 

In this section we use the entropy curves derived in the pre- 
vious section for the unit filled lattice to calculate how the 
temperature of a dilute 3D Bose-Einstein condensate (BEC) 
changes as it is adiabatically loaded into a deep optical lattice. 
Ideally the adiabatic loading process will transfer a T = 
BEC into a perfect Mott Insulator (MI), however condensates 
can not be created at = and it is important to know the re- 
lation between final and initial temperatures. Calculations for 
an ideal bosonic gas 0] demonstrate that for typical tempera- 
tures at which a BEC is created in the laboratory, adiabatically 
ramping up the lattice has the desirable effect of cooling the 
system. On the other hand, drastic changes in the energy spec- 
trum (the opening up of a gap) induced by interactions modify 
this ideal situation jj] and in the interacting case atoms can be 
instead heated during the loading. 

In order to calculate the change in the temperature due to 
the loading, we first calculate the entropy as a function of 
temperature of a dilute uniform BEC of ^^Rb atoms by using 
Bogoliubov theory. The Bogoliubov approximation is good 
for a dilute gas as it assumes that quantum fluctuations intro- 
duced by interactions are small and treats them as a perturba- 
tion. The quartic term in the interacting many-body hamilto- 
nian is approximated by a quadratic one which can be exactly 
diagonalized fl3lfl4ll . This procedure yields a quasi-particle 

excitation spectrum given by Cp = \J i^fY + 2unejj. Here 

Ep — p^/2m are single particle energies, u ~ Anh^as/m, m 
is the atomic mass and n is the gas density . 

Using this quasi-particle spectrum in the Bose distribution 
function of the excited states, /(cp) = [e'^'^p — 1]^^, one can 
evaluate the entropy of the gas given by 

^k.,„,.=o = k 5^{/?6p/(ep) - ln[l - e^^^]}. (17) 
p 

Using Eq.dlVli we numerically calculate the entropy of the sys- 
tem for a given initial temperature T^. Assuming the entropy 
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FIG. 2: (color online) Top: Tf vs Ti (in units of Er) for differ- 
ent final lattice depths Vf. Here, we assume adiabatic loading in 
the limit J = 0. The dashed(red), dot-dashed(blue), and long- 
dashed(black) lines are for Vf — 10, 20 and 30Er respectively. The 
continuous(grey) lines are calculated for the various lattice depths 
from Eq.(19). The dotted line is the identity, Tf = T. Bottom: Av- 
erage density of unit filled cells ni as a function of T (in units of 

Er). 



during the adiabatic process is kept constant, to evaluate Tf 
for a given Ti we solve the equation, 

5(r.)|K.,„,,.0 = 5(T/)|y,,„,,=y,. (18) 

We evaluate the right hand side of this equality assuming that 
the final lattice depth, Vf, is large enough that we can ne- 
glect terms proportional to J in the Hamiltonian. We use 
the expression for the entropy derived in the previous section, 
Eq. (12), together with Eqs. (14) and (15). 

The results of these calculations are shown in Fig |2] where 
we plot Tf vs Ti for three different final lattice depths, 
Vj: /Er = 10 (dashed Une), 20 (dot-dashed line) and 30 (long- 
dashed line). In the plot both Tf and Ti are given in recoil 
units Eji. As a reference, the critical BEC temperature for 
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an ideal bosonic gas (which for a a dilute gas is only slightly 
affected by interactions) in recoil units is kT^ « Q.67Ej^. 

For kTi > 0.05£'/? the final temperature scales linearly 
withT,;: 



fcr, 



/ 



u 

3^ 



(19) 



In Fig. 121 Eq.(I19> is plotted with a gray line for the various 
final lattice depths. 

In contrast to the non-interacting case, where for kTi < 
O.bEu the system is always cooled when loading into the 
lattice 01, here interactions can heat the atomic sample for 
low enough initial temperatures. For reference in Fig. 2, we 
show the line T/ T^. One finds a temperature r'*''"*(F/), 
(determined from the intersection of the Tf = Ti line with 
the other curves) below which the system heats upon load- 
ing into a lattice of depth V/. From the linear approxima- 
tion one finds that T'"'"* increases with U as kT'^''°-\Vf) sa 
0.177?7(3 — U /Er)^^. Because U scales as a power law with 
the lattice depth a larger Vf implies a larger r'"^°*(y^) 
and so a larger heating region. Note that for the shallowest 
lattice in consideration, Vf/En = 10, kT^""'* < 0.05 and 
therefore the linear approximation does not estimate it accu- 
rately. Fig. 2 also shows a very rapid increase in the tem- 
perature close to Ti = 0. This drastic increase is due to the 
low temperature plateau induced by the gap that opens in the 
insulating phase. 

To quantify the particle-hole excitations and give an idea of 
how far from the target ground state the system is after the 
loading process, we also plot TTi vs Ti in the bottom panel of 
Fig. 2 . In the plot, ni is calculated from Eq.(16). We found 
that to a very good approximation 



ni(Ti] 



1 — exp 



2kT/ER + 0.354, 



(20) 



Note that in the J — hmit, ni depends exclusively on f3U 
and thus as long as the final lattice depth is large enough to 
make the J = approximation valid, ni is independent of 
the final lattice depth. The exponential suppression of multi- 
ple occupied states in the entropy plateau explains why even 
though the final temperature increases rapidly near Tj = 0, 
this is not reflected as a rapid decrease of rii. For the largest 
initial temperature displayed in the plot, kTi/En w /2, 
the final temperature reached in units of U is kTf /U w 0.17 
and Til ~ 0.9. Thus, the fidelity of the target state has been 
degraded to less than 90%. In Fig.l one also observes that 
ni « 0.9 at kTf/U « 0.17 and that most of the loss of fi- 
delity is due to particle-hole excitations as nr>3 ~ 0. 



C. Finite size effects 

In recent experiments by loading a BEC into a tight two- 
dimensional optical lattice, an a rray of quasi-one dimensional 
tubes has been created ll El [iH [13. Ill- The number of 
atoms in each tube is of the order of less than 10^ and there- 
fore the assumption of being in the thermodynamic Umit is no 
longer valid for these systems. 



The thermodynamic limit assumption used in the previous 
section allowed us to derived thermodynamic properties with- 
out restricting the Hilbert space in consideration. Thus, within 
the one-band approximation, these expressions were valid for 
any temperature. However, if the size of the system is finite, 
number fluctuations AiV must be included and to derive ex- 
pressions valid for arbitrary temperatures could be difficult. In 
this section we calculate finite size corrections by restricting 
the temperature to kT < U /2. At such temperatures Fig. [l] 
shows that only states with at most two atoms per site are rele- 
vant so one can restrict the Hilbert space to include only states 
with at most two atoms per site. 

Setting nr>2 = and M = in Eq.(|3, the partition 
function (at zero order in J ) can be explicitly written as: 



NI 



L7V/2J 



= e-P^/^cos{'.N)&^''\\ef^^l% (21) 

where Ci™' [x] are Gegenbauer polynomials j^l- 

In Fig |3l (left panel) we study the effect of finite atom num- 
ber on the entropy. We show the entropy per particle as a 
function of temperature for systems with = 50 (green 
dotted fine), N = 100 (blue dashed fine) and N = 1000 
(red dash-dotted line). For comparison purposes we also plot 
with a (black) solid line the entropy calculated using Eqs. J14t 
and M5\ . which were derived in the thermodynamic limit. It 
can be observed that for N = 1000 the thermodynamic limit 
is almost reached (nearly indistinguishable from the thermo- 
dynamic limit). Finite size effects decrease the entropy per 
particle and thus tend to increase the final temperature during 
the adiabatic loading. 

Furthermore, in the right panel we also compare Eq. Mil 
with the entropy calculated by restricting the Hilbert space 
even more and including only one-particle-hole (1-ph). 1- 
ph excitations are the lowest lying excitations which corre- 
spond to states that have one site with two atoms, one with 
zero atoms and one atom in every other site, i.e. {nr}u = 
{1, TV - 2, 1, 0, . . . , 0}. There are N{N - 1) different particle 
hole excitations all with energy U . If the entropy is calculated 
taking into account only 1-ph excitations one gets an expres- 
sion to zeroth order in J given by: 



S _ \n[l + N{N -l)e-'^^] 
k 



N 



(3U{N -l)e-'^^ 
l + N{N - l)e-PU- 



(22) 



The right panel shows that as long as the temperature is 
below kT ^ O.IC/ and the number of wells is of order 10^ or 
less, Ea.( l22t gives a very good approximation for the entropy 
per particle. 



D. Finite J corrections 

In the previous section for simplicity we worked out the 
thermodynamic quantities assuming J = 0. However, if the 
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FIG. 3: (color online) Left: Entropy per particle 5 as a function of the temperature T (in units of U) for a unit filled lattice in the J = limit 
and different number of atoms A'^. The solid line shows S calculated in the thermodynamic limit using Eq. <14> and <15> . The dash-dotted(red), 
dashed(blue), and dotted(green) lines correspond to A'^ = 1000, 100 and 50, respectively. For these curves, S is restricted to the ph subspace 
(see Eq. <21> '). Right: S vs T (in units of U) for A'^ = 100. The dashed and solid lines are the entropy calculated in the ph and l-ph(see Eq I22> 
subspaces, respectively. 



final lattice is not deep enough, finite J corrections should be 
taken into account. In this section we study how these cor- 
rections can help to cool the unit filled lattice during adiabatic 
loading. 

In the J = limit all thermodynamic quantities are inde- 
pendent of the dimensionality of the system. On the other 
hand, for finite J the dimensionality becomes important. In- 
cluding J in the problem largely complicates the calculations 
as number Fock states are no longer eigenstates of the many- 
body Hamiltonian and many degeneracies are lifted. For sim- 
plicity, in our calculations we will focus on the ID case and 
assume periodic boundary conditions. We will also limit our 
calculations to systems with less than 10^ atoms and tempera- 
tures low enough (kT <C O.IU) so it is possible to restrict the 
Hilbert space to include only 1-ph excitations. 

To find first order corrections to the N{N 1) low lying 
excited states we must diagonalize the kinetic energy Hamil- 
tonian within the 1-ph subspace. For 1-D systems this diag- 
onalization yields the following approximated expression for 
the eigenenergiesQ^ 

i.«=t/-4Jcos(^)cos(f), (23) 

Where r = 1, ... - 1 and i? = 0, ... TV - 1. Using these 
eigenenergies to evaluate the entropy per particle one obtains 
the following expression: 

^ « ^""^grpf.r [/g(2J/?)-f/o(2J/3)Ji(2J/3)] ^_,^ 



with 

Z = 1 + N{N -l)e-'^^Io{2J(3), (24) 
where /„(a;) are modified Bessel functions of the first kind 

To derive Eq.( I23> . we assumed similar effective tunneling 
energies for the extra particle and the hole. This is not exact, 
especially for a unit filled lattice, g = 1, since the effective 
hopping energy for the particles and holes goes like J{g + 1) 
and Jg respectively. However, we find by comparisons with 
the exact diagonalization of the Hamiltonian that for observ- 
ables such as the partition function which involves summing 
over all the 1— ph excitations, this assumption compensates 
higher order corrections in J/U neglected to first order It 
even gives a better expression for the entropy of the many- 
body system than the one calculated by using the spectrum 
obtained by exact diagonalization in the 1— ph subspace. 

We show the validity of Eq.(|3^ in Figl^where we compare 
its predictions ( plotted with solid lines) with the entropy cal- 
culated by exact diagonalization of the Bose-Hubbard Hamil- 
tonian for different values of J/U assuming a system with 
= 10 atoms. In the plot we use a dot-dashed line for 
J/U = 0.1 and a dashed line for J/U = 0.01 . Even for 
the case J/U = 0.1 we see the analytic solution reproduces 
very well the exact solution, especially at low temperatures. 
AtfcT > 0. lie/ higher order corrections is J/C/ become more 
important. 

We now use Eq.(I24> to study larger systems where an ex- 
act diagonalization is not possible. Even though we expect 
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0.025 0.05 0.075 0.1 0.125 0.15 
kT/U 

FIG. 4: (color online) Finite J corrections: The dash-dotted and bro- 
ken lines correspond to the entropy per particle vs T (in units of U) 
calculated by numerical diagonalizations of the Hamiltonian for sys- 
tems with A'^ = 10 atoms and J/U=0.l and 0.01 respectively. The 
corresponding solid lines show the entropy per particle calculated 
from Ea.i24l. 



finite J corrections to become important at lower tempera- 
tures for larger systems, we consider that for systems with 
less than 10^ atoms, small values of J/U and within the low 
temperature restriction, Eq.( I24> can still give a fair descrip- 
tion of the entropy. In Fig|5lwe show the effect of finite J 
corrections on the final temperature of a system of 100 *^Rb 
atoms when adiabatically loaded. For the calculations, we 
fix the transverse lattice confinement to Vx = Vy = 30Efj, 
assume d = 405 nm and vary the axial lattice depth. We 
show the cases Vz = ^Er with a yellow dash-dot-dotted line, 
Vz = 7 Ell with a red dashed line, Vz = 9Er with blue dot- 
ted line and Vz = llEji with a green dash-dotted line. For 
these lattice depths, the single-band approximation is always 
valid and the energies J and U both vary so that their ratio 
decreases as J/U = {0.12, 0.07, 0.04, 0.02} respectively. We 
also plot for comparisons purposes the J/U = case with a 
solid black line. 

Fig|5lshows that finite J corrections decrease the final avail- 
able temperature of the sample. These corrections are impor- 
tant for shallower lattices, as they decrease the final tempera- 
ture with respect to the J = case by about 30%. For lattices 
deeper than Vz = lli^i? the corrections are very small. 

The decrease in the final temperature induced by J can 
be qualitatively understood in terms of the modifications that 
hopping makes to the eigenenergies of the system. J breaks 
the degeneracy in the 1— ph, leading to a quasi-band whose 
width is proportional to J. As J increases the energy of the 
lowest excited state decreases accordingly, while the ground 
state is only shifted by an amount proportional to j"^ /U . The 
lowest energy excitations then lie closer to the ground state 
and become accessible at lower temperatures. As a conse- 
quence, the the entropy increases (and thus Tf decreases) with 



respect to the J = case. 

Following the same lines of reasoning the entropy should 
exhibit a maximum at the critical point associated with the 
Mott insulator transition, since at this point an avoided cross- 
ing takes place. We confirmed this intuitive idea with exact 
numerical diagonalization of small systems. For the trans- 
lationally invariant case, we expect the entropy to become 
sharply peaked at the transition with increasing N and this 
could be an important limitation for adiabatically loading 
atoms. However, as we will discuss later, the external har- 
monic confinement present in most experiment prevents a 
sharp Mott insulator transition and can help to decrease the 
adiabaticity loading time within accessible experimental time 
scales. 

In this section we focused on the effect of finite J correc- 
tions in ID systems. For higher dimensions, we expect that 
finite J corrections help to cool the system even more, since 
the effective tunneling rate that enters in the entropy scales 
with the number of nearest neighbors and thus becomes larger 
for higher dimensions. 



IV. HARMONIC CONFINEMENT: Vi = Qi^ 

For simplicity in our analysis we consider a ID system 
which can be studied using standard fermionization tech- 
niques iPqI. These techniques allow us to map the com- 
plex strongly correlated bosonic gas into a non-interacting 
fermionic one. We choose our parameters so that they closely 
resemble the experimental ones used in Ref.|01. Specifically 
we use transverse lattice depths of Vx = Vy = 27 Eb, created 
by lasers with wavelength = 823 nm and an axial lat- 
tice depth of Vz = I'&.hEji created by a laser of wavelength 
\z ~ 854 nm. We set the axial frequency of the ID gas to 
LOz ~ In y. 60 Hz and the number of atoms to iV = 19 (this 
was the mean number of atoms in the central tube of the ex- 
periment). For these parameters the ratio U / J ^ 205 and 
Vl/J = 0.28 with f7 = l/8mcj2A2. The ground state of the 
system corresponds to a MI with iV unit filled sites at the trap 
center (see FigQbottom panel). We compare the thermody- 
namic properties of this system with the ones of a translation- 
ally invariant system in the MI state, with the same number of 
atoms (A^ ~ M ~ 19) and same ratio U /J. As we described 
in the previous section, for homogeneous systems the finite 
,/ corrections for the J/U ^ 0.005 ratio in consideration are 
very small, and for temperatures below kT/U < 0.1 they can 
be neglected. On the contrary, when the parabolic confine- 
ment is present, taking into account the kinetic and trapping 
energy corrections is crucial for a proper description of the 
low temperature properties of the system. Unlike the spatially 
invariant case, where the lowest lying excitations in the MI 
phase are 1-ph excitations which have an energy cost of order 
U , in the trapped case, within the parameter regime in consid- 
eration, there are always lower-lying excitations induced by 
atoms tunneling out from the central core and leaving holes 
inside it. We refer to these excitations as n-hole (n-h) exci- 
tations. These "surface" n-h excitations must be included in 
the trapped system because of the reservoir of empty sites sur- 
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FIG. 5: (color online) Tf (in units of U) vs Ti (in units of Er) curves calculated using Ea. <24> for a system with A'^ = 100 atoms and different 
values values of J/U: dash-dot-dotted(yellow) line J/U = 0.12(V3 = 5Er), dashed(red) line J/U = 0.07(V3 = 7Er), dotted(blue) line 
J/U = 0.04(14 ~ and dash-dotted (green) line J/U — 0.02(14 = IIEr). The solid(black) lines is shown for comparison purposes 

and corresponds to the limiting case J — 0. 



rounding the central core, which introduces an extra source of 
delocalization. For example the lowest lying hole excitations 
correspond to the 1-h excitations created when a hole tunnels 
into one of the most externally occupied sites. They have en- 
ergy cost QN, which for the parameters in Ref. Q] is 40 times 
smaller than U. 

For a system in arbitrary dimensions, it is complicated to 
properly include n-h excitations in the calculations of ther- 
modynamic properties. For ID systems, however, the Bose- 
Fermi mapping allows us to include them in a very simple 
way. Nevertheless, because fermionization techniques neglect 
multi-occupied wells in the system we have to restrict the 
analysis to temperatures at which no multiple occupied states 
are populated (kT < O.IU, see also Ref.01). The results are 
shown in Fig|6l where we plot the final temperature of the 
sample as a function of a given initial entropy S. In the plot 
we also show the results for the corresponding translationally 
invariant system. 

The most important observation is that instead of the sud- 
den temperature increase at = (or flat S vs T plateau 
induced by the gap), in the trapped case the temperature in- 
creases slowly and almost hnearly with S: 

^ « A (^) , (25) 
A = 5fc(^)', (26) 



with Tp the Fermi temperature. The linear behavior is char- 
acteristic of low temperature degenerate Fermi gases and the 
proportionality constant A depends on the density of states of 
the system. For this particular case, A can be estimated as- 
suming a box-like dispersion spectrum, £'„ = rin^. For the 
parameter regime in consideration this assumption is valid for 
the modes close to the Fermi energy, which are the relevant 
ones at low temperature (Ref. l22ll ). Using this box-like spec- 
trum it is possible to show that for ft < kT ^ kTp (where 
the first assumption allows a semiclassical approximation) A 
is given by Ea.(l26>. In Fig|6lthe blue-dotted line corresponds 
to this linear solution and it can be seen that it gives a fair de- 
scription of the entropy in the low temperature regime. It is 
interesting to point out that the slower increase in entropy as 
a function of temperature in the homogeneous system com- 
pared to the trapped one is a particular effect induced by in- 
teractions. In the non-interacting case the opposite behavior 
is observed; for an homogeneous system S*'' cx (T/Tc)^^^ 
whereas for the trapped system 5" cx (T/Tc)^ so if T < Tc 
then S'^ > S'^. Here D is the dimensionality of the system 
and Tc the critical condensation temperature. 

In a typical experiment the sample is prepared by forming 
a BEC in a magnetic trap. Therefore a good estimate of the 
initial entropy is given by IIH : 

S^k f 4C(4)/C(3)t3 + ,^1^2 (1 _ t3)2/5\ (27) 
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FIG. 6: (color online) Final temperature (in units of U) vs initial en- 
tropy per particle. The solid(black) line corresponds to the trapped 
system with the parameters chosen to closely reproduce the experi- 
mental set-up of Refj^l- The dotted(blue) line is the analytic solu- 
tion for a system of fermions at low temperature assuming a box-like 
spectrum and the dashed(red) line is the entropy for the correspon- 
dent homogeneous system calculated from Eq. <22> 

where 77 = a(7Vi/6as/afto)2/5, with a = 12 « 

1.57 and a/,,, = ^hj (muj) the mean harmonic oscillator 
length (ZU = :^ujxijjyijj^. The parameter r\ takes into account 
the main corrections to the entropy due to interactions. In 
the above equation t = T/T^ with = T°{1 - 0.4377^/2) 
the critical temperature for condensation and the critical 
temperature for the ideal trapped gas kT^ = fiIU{N/({3)y^^. 
The term proportional to 77^/^ accounts for the small shift in 
the critical temperature induced by interactions. For typical 
experimental parameters 77 ranges from 0.35 to 0.4. 

If one assumes a TV = 10^'^ atoms, aJ/(27r) = 60 - 120 
Hz and a very cool initial sample, t ^ 0.2, one obtains that in 
typical experiments the initial entropy per particle of the sys- 
tem is not smaller than S/k > 0.1. Fig|6lshows then that the 
reduction of the final temperature during the adiabatic loading 
induced by the trap can be significant. In turn, this suggests 
that the presence of the magnetic confinement is going to be 
crucial in the practical realization of schemes for lattice-based 
quantum computation. 

To emphasize this point, in Fig. (top panel) we plot the 
mean occupancy of some lattice sites as a function of the ini- 
tial entropy per particle. It should be noted that for the number 
of atoms in consideration the edge of the cloud at T = is at 
j = (A^— l)/2 = 9. For comparison purposes we also plot ni 
calculated for the correspondent spatially homogeneous sys- 
tem. The plots shows that for the central lattice sites there is 
almost 100% fidelity to have one atom per site for the range 
of initial entropies in consideration. Fluctuations are only im- 
portant at the edge of the cloud and if one excludes these ex- 
tremal sites the fidelity in the trap case remains always higher 
than the fidelity in absence of the trap. In the bottom panel we 
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FIG. 7: (color online) Top: Local on-site probability of having 
unit filling ni as a function of the entropy per particle S, and for 
a few values of the site-index j. The dashed(red), dotted(blue), dot- 
dashed(black), and long-dashed(black) lines correspond to the sites 
with j — 0,4,7 and 8, respectively. For comparison purposes, we 
also plot ni calculated for the correspondent spatially homogeneous 
system (solid line). Bottom: Density profile for the trapped system 
in consideration at S/k — (dashed line) and S/k = 0.25 (solid 
line). 



also show the density profile for S* = and compare it with 
the one at S/k = 0.25. It is clear in the plot that the central 
lattice sites remain almost with one atom per site. 

The considerations made here are for adiabatic changes to 
the lattice, and therefore represent a lower bound on the final 
temperature, assuming the entropy is fixed. How quickly the 
lattice can be changed and remain adiabatic is a separate is- 
sue, but we point out that for systems with finite number of 
atoms confined by an external trap there is not a sharp super- 
fluid/insulator phase transition, which should relax the adia- 
baticity requirements when passing through the transition re- 
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gion. A proper adjustment of the harmonic confinement dur- 
ing the loading process could reduce the time scales required 
for adiabaticity to be experimentally realizable. 

V. CONCLUSIONS 

In this paper we calculated entropy-temperature curves for 
interacting bosons in unit filled optical lattices and we used 
them to understand how adiabatic changes in the lattice depth 
affect the temperature of the system. 

For the uniform system, we have derived analytic expres- 
sions for the thermodynamic quantities in the J /U = case 
and we used them to identify the regimes wherein adiabati- 
cally changing the lattice depth will cause heating or cooling 
of the atomic sample in the case of a unit filled lattice. We 
have shown that the heating is mainly induced by the gapped 
excitation spectrum characteristic of the insulator phase. By 
considering finite size effects and finite J corrections we have 
shown that the former leads to increased the heating of the 
atoms, the latter tend to reduce it. 

Finally, we have discussed the spatially inhomogeneous 
system confined in a parabolic potential and we have shown 



that the presence of the trap reduces significantly the final 
available temperature of the atoms due to the low-energy sur- 
face excitations always present in trapped systems. 

The fact that the harmonic confinement turns out to be 
clearly a desirable experimental tool for reducing temperature 
in the lattice is an important finding which should be taken 
into account in the ongoing experimental and theoretical ef- 
forts aimed at using the Mott Insulator transition as a means 
to initialize a register for neutral atom quantum computation. 
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